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, Abstract 

! An approximate formula connecting the true and the HBT ho- 

^> I mogeneity regions in multparticle production processes is derived. It 

■ impUes that when calculating the HBT radii one should use the cen- 
. ter of mass systems of the pairs rather than the now popular LCMS 
I system. A discussion of several simple examples clarifies the poten- 

■ tial and limitations of the HBT method. The even cumulants of the 
. X-distribution, including the HBT radii, can be determined for each 
I homogeneity region, but the relative positions of the homogeneity re- 
' gions are unconstrained. This makes the HBT radii of little use for 

calculating quantities dependent on the interparticle interactions in 
coordinate space. 



X 



PACS numbers 25.75. Gz, 13.65.+i 

Bose-Einstein correlations, interaction region determination. 



Almost fifty years ago, in a famous paper known as GGLP [Ij, a method 
of using momentum distributions for pairs of identical pions to estimate the 
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sizes of the interaction regions, i.e. of the regions where the hadrons are 
produced in multiparticle production processes, has been described. Such 
radii, determined from momentum distributions of identical particles, have 
been later called, not very appropriately j2j, HBT radii. GGLP assumed 
that hadron production happens instantly and simultaneously at some time 
t = and that there are no correlations between the momenta of the hadrons 
and their production points. Considering simultaneously the momentum of 
a particle and its production point implies a quasiclassical approximation, 
but it can be made plausible [3J that this approximation is good for heavy 
ion scattering and acceptable also for the other multiparticle production pro- 
cesses. Let us denote by pi and p2 the four-momenta of the two identical 
pions in the paii0, by Xi and X2 the space-time positions of their production 
points and introduce the notation 

K = ^{pi+P2); q=Pi-P2] X = ]^{xi+ X2); y = Xi-X2. (1) 
Note that Pi = P2 = ""^^ implies 



Ko = -^^El + + ^iEl + iq2)2 - (K ■ q)^; go = /3 ■ q, (2) 

where Ek = \/K^~Tm^ equals -ft^o at q = 0, and (3 = ^ is the velocity of 
the pair. GGLP found Rhbt from the distribution of g^. 

Some twenty five year later Pratt ^ described a model with an exploding 
source. In this model, the momentum distribution of particles depends on 
the production point. Due to this correlation, the HBT radius differs signifi- 
cantly from the true radius, which is known from the input. For fixed K the 
result was that, while the true radius R does not depend on |K|, the HBT ra- 
dius decreases from y^-R for K = to zero when |K| tends to infinity. This 

finding was generalized by Bowler [5], who pointed out that in general, when- 
ever there are strong correlations between momenta and production points, 
Rhbt measures only the region where particles of similar momentum are 
produced. Sinyukov [6], [7] associated these regions with the homogeneity 

^We discuss pions for definiteness, but the discussion applies to any spin zero bosons. 
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regions considered in hydrodynamics and the name homogeneity region got 
generally accepted. 

Most models in their simplest form, i.e. without final state interactions, 
resonance decays etc., find the ra-particle momentum distributions for iden- 
tical pions from formulae equivalent to [8j 

n 

^(Pl, •••,?«) =Cn^Wp{Pj-.VQj). (3) 

Q i=i 

where the summation is over all the permutations j — > Qj of the indices 
j = l,...,n, Cn are normalization constants and p(Pi;P2) is some time 
independent, single particle density matrix. Thus p(p]^;p2) determines all 
the momentum distributions. The best [9J quantum mechanical analogue of 
the classical phase space density is the Wigner function iy(X, K), which is 
a Fourier transform of p(Pi; P2)- Following most models, we assume further 
that, W(K, K) defines the HBT X- distribution for the homogeneity region 
corresponding to a given K. 

Measuring all the momentum distributions is not enough to determine 
p(PijP2)- As easily seen from ([2]), the observable momentum distributions 
do not change [10], [H], under the transformatioij^ 

P(pi; P2) - P/(Pi; P2) = e^[^(^^)-^(^^)ip(Pi; P2), (4) 

where f{p) is an arbitrary real- valued function of the four- vector p, and 
Po = \/p'^ + m?- This transformation does not affect the HBT radio of the 
homogeneity regions, but it can shift and deform these regions. 

An important concept introduced in [1] is the emission function S{X^ K). 
Assuming chaotic sources, i.e. no interference between particles produced 
at different moments of time, this is related to the density matrix in the 
momentum representation by the formula^ 

^The full group of transformations includes also p(pi; P2) P*(Pi; P2)i but this corre- 
sponds to the space inversion of the interaction region and is of no interest in the present 
context. 

■^More generally, all the even cumulants of M^(X,K), at K fixed, remain unchanged. 
The proof is a simple extension of the proof given for the HBT radii in [11] . 

^Sometimes formulae equivalent to ([5]) with p(pi;p2) replaced by p*(pi;p2) are used. 
The advantage of the present convention is that for S{X, K) = 6{t)g{'K, K) function g is 
just the Wigner function. 
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p(Pi;P2)= d^XS{X,K)e 



MX 



(5) 



In the spirit of the quasiclassical approximation, using the well tested (cf. 
e.g. [12] and references given there) mass shell approximation 



one can interpret S{X, K) as the time-dependent distribution of the pairs of 
vectors {X, K}. The time independent homogeneity region for each K can 
be obtained by integrating over time 



where the notation stresses that we are interested in the distribution of X 
for given K. At this point the mass shell approximation ([6]) is necessary in 
order to make Kq, and consequently S{X;K) and p(X|K), independent of 

q 

Further we will call p(X|K) true distribution in order to distinguish it 
from the HBT distribution given by the Wigner function. Invariance (jll) 
means that for a given set of momentum distributions there is a variety of 
HBT X-distributions which all correspond to the same fit to the data. One 
way of choosing among them is to specify the emission function. 

It was soon noticed [13] that for a given p(p^;p2) there is an infinity 
of solutions for S and that, in particular, a small long-lived source may be 
undistinguishable from a large short-lived source. Thus, the information 
about interaction regions obtained from momentum measurements is rather 
incomplete. 

The standard relation between the density matrix and the Wigner func- 
tion yields from ^ 



where go has been eliminated using ([2]). We invoke now the mass shell ap- 
proximation ([6]), which makes S independent of q. Thus, the integration d^q 
gives (27r)^(5(X — X' -|- f3t) and the integration d^X' can be performed. The 
result is 




(6) 




(7) 
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Wi'K.K) = Jdt S{:S. + /3t,t,K). (9) 

The comparison of the true interaction region with the HBT one reduces to 
the comparison of the integrals in ([7j) and (Q. 

The integrals in ([7]) and (j9), and consequently the true and HBT homo- 
geneity regions, coincide onljQ for /3 = 0. This can be achieved by working 
with pairs which have the same velocity f3 and using the reference frame 
where f3 = 0. 

We summarize our findings: 

• The HBT results are credible only for single homogeneity regions, i.e. 
at given K. 

• For each homogeneity region, only the pairs with the same velocity, 
i.e. with the same K, are considered and one should use the reference 
frame, where this velocity is zero. There, the true and the HBT homo- 
geneity regions coincide which is not the case for other frames, as for 
instance for the now popular (cf e.g. the review [H]) LCMS frame. An 
additional pragmatic argument in favor of the rest frame is that the 
effects of the final state interactions are the simplest there (cf. e.g. [17] 
and the references given there). 

• The true radii of the homogeneity regions, more generally all the even 
cumulants of their X-distributions, can be obtained from the HBT 
analysis of the momentum distributions. For instance, when p(X|K) = 
p(— X|K) all the odd cumulants vanish and, therefore, p(X|K) can be 
measured by the HBT method. 

• The space distribution of the centers of the homogeneity regions (X)(K) 
is unconstrained by the HBT analysis, but once it is fixed, all the func- 
tions p(X|K) can be determined [10], [TT] . 

Let us discuss some simple examples illustrating these features of the 
HBT method. In order to avoid unessential complications, we will consider 
one space dimension and very simple emission functions constructed from 
5-functions and step functions. They violate the Heisenberg uncertainty re- 
lations, but all these calculations can be repeated in three dimensions using 
Gaussians and the results are qualitatively the same. 



^Except for the trivial case when S(X, K) ^ 6{t). 
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Let us first consider the emission function 

K- to, a, 6, k) = 6{t - to)e,,6(X)5(K - k), (10) 

where to, a < b and k, are real constants. Function Qa,b(yX) = ^ for a < x < 6 
and zero outside this interval. For this emission function 

p(p,;p,) = ^(K-«:)e--(^-/5*«)^|^^ (11) 
1^(X,K) = e„,b(X + /3to)<5(K-/€). (12) 

It is seen that the segment a < X < 6 got shifted to a — f3to < X < 6 — f3to- 
Acording to (ITT]) this shift is due to the phase factor in the density matrix. 
Since in the formula for the two-particle correlation function only the absolute 
value of p appears, experiment is blind to such shifts. This remains true also 
when many-particle correlation functions are measured, even when all the 
measurements are performed with perfect precision [10], [11]. The fact that 
the position of the center of the interaction region cannot be found from the 
HBT analysis of the momentum distributions is, of course, well known (see 
e.g. [12]), but one should also keep in mind that when the emission function 
is given the density matrix, the Wigner function and, consequently, all the 
HBT homogeneity regions are unambiguously defined. 
As the next example we take 

S^iX, K) = ^ [S,iX, K; to, -a, 0, /ti) + S,{X, K; to, 0, a, f^^)] (13) 

Using the previous example, it is seen that the segment —a < X < a is 
mapped onto two segments: —a — fi^to < X < —fB^to and — /32^o < X < 
a — /32^o- Let us choose 

f{p) = -b ■ p - icp^ (14) 

where b is an arbitrary vector and c an arbitrary constant, and make the 
transformation (j4]). Each HBT homogeneity region gets shifted by b -|- cK. 
For Ki 7^ K2, by a suitable choice of c, one can obtain any prescribed distance 
between the centers of the two segments. The true length of the interaction 
region, as seen from S, is 2a. The length of the transformed interaction re- 
gion, which follows just as well from the data on momentum distributions. 
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can be any number not smaller than a. What is the way out? One has to 
invoke the homogeneity regions. For each K separately, the length of the 
segment where particles with this value of K are produced, i.e. of the homo- 
geneity region, is reproduced correctly. The positioning of the homogeneity 
regions corresponding to different values of K, however, is beyond control 
when the momentum distributions are the only input. In three dimensions it 
is easy to prove [TU], [H] that by a suitable choice of function f{p) in (jl]), the 
positions of the centers of the homogeneity regions (X)(K) can be changed 
into (X)(K) + g(K), where g is an arbitrary differentiable function of K. 

As an amusing example in three dimensions let us consider the models 
where ^ AX^, and A is a constant [15], [TB]. They correspond to emission 
functions 

S{X,K) = 6%X- f3t)S{X,t,K), (15) 

where S is some function which does not affect the singularity introduced by 
the 6^. The corresponding Wigner function is 

W(K,K) = 5%X.) jdtS{f3t,t,K). (16) 

Thus, the HBT interaction region reduces to one point. To be sure: such 
models, when properly used, are quite successful, but finding the interac- 
tion region from the density matrix in the momentum representation is their 
misuse. 

As our last model consider the emission function 

S{X, K) = eo,a(t)5(X)5(K - n). (17) 
This corresponds to the Wigner function 

M^(X, K) = ieo,„(-X//3)5(K - k) = 0_^^^q(X)5(K - k). (18) 

The length of the true interaction region is zero, while the length of the HBT 
interaction region is /3a. Thus, in order to get the true length from the HBT 
analysis one must use the reference frame where K = 0. We conclude that, 
working with all the pairs which have a given velocity, one should measure 
their homogeneity region in their rest frame. Let us consider some further 
implications of our analysis: 
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If there are no position-momentum correlations, all the homogeneity re- 
gions, measured in the respective rest frames, are equivalent to each other 
and to the overall interaction region. The ambiguity (jlj) reduces to a lack 
of information about the position of the center of the interaction region [11] . 
Thus, in this case the HBT method works very well. One should keep in 
mind however that, as illustrated by our second model, the independence of 
the homogeneity region on K is a necessary, but not a sufficient condition 
for the absence of position-momentum correlations. 

If there are no interparticle interactions in coordinate space, the HBT 
homogeneity regions may be all we need. E.g. the entropy of a gas of 
noninteracting particles (in the quasiclassical approximation) depends on the 
accessible phase space volume. This can be calculated by integrating over 
coordinate space at fixed momentum and then integrating over momenta. 
The first integration gives just the HBT volume of the homogeneity region, 
though one should keep in mind the ambiguity (jl]). If there are interparticle 
interactions in coordinate space, however, it may make a lot of difference 
whether the homogeneity regions are on top of each other, or scattered over 
space. Thus, the use of the HBT radii to calculate the entropy of a gas of 
interacting particles, or of their mean free paths, is risky. 

For f3 = the true space density Jdt S{X, K) can be unambiguously 
obtained from p(Pi, P2) by inverse Fourier transformation. Thus, all the am- 
biguity in the determination of the homogeneity regions, including their abso- 
lute positions, results from and disappears when f{p) is fixed. Conversely, 
when the true space density is known, it yields unambiguously p(pi; P2) and 
thus it fixes f{p). 

The ambiguity in the determination of the space-time density S follows 
from the irreversibility of transformation ^ and persists even when f{p) is 
fixed, i.e. when p(p]^;p2) is known. 
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